AD-A171  822  RANDOM  CIRCLES  AND  FIELDS  ON  CIRCLES(U)  NORTHMESTERN  1/1 

UNIV  EVANSTON  IL  E  CINLAR  ET  AL.  AAV  88 
AF0SR-TR-8S-8S57 


UNCLASSIFIED 


F/Q  12/1 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la  W E P^RT  SECU m TV  Ct ASS4 R C A T ION 


2a  security  classification  authority 

M—— - 

2b  DEC  LASS  IF  1C  ATI  ON/D  OWNGRAQlNG  SCHE 


NG  ORGANIZATION  REPORT  NUMBER(S) 


AD-A171  833 


3 1  LI T  Y  OF  REPORT 


pproved  for  public  release 
■  unclassified  distribution  unlimited. 


6.  NAME  OF  PERFORMING  ORGANIZATION 

Northwestern  University 


6c.  ADDRESS  (City.  Slate  and  ZIP  Code I 

Evanston ,  IL  60201 


b.  OFFICE  SYMBOL  7*.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable  i 

Air  Force  Office  of  Scientific  Research 


7b.  ADORESS  (City.  State  and  ZIP  Code ) 


Washington,  D.C. 


.  NAME  OF  FUNOING/SPONSOHING 
ORGANIZATION 


4S£ 


18b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

(If  applicable!  1 

|  AFOSR  Grant  No.  82-0189 


10.  SOURCE  OF  FUNOING  NOS. 


PROGRAM 
ELEMENT  NO. 


11.  TITLE  ( Include  Security  Classification ) 

Random  Circles  and  Fields  on  Circles 


12.  PERSONAL  AUTHOR(S) 

Erhan  QINLAR  and  J.  G.  WANG 


13b.  TIME  COVERED 


PROJECT 

NO. 

TASK 

NO 

WORK  UNIT 
NO 


17. 


FIELD 


COSATI  COOES 


18-  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number > 

Gaussian  fields,  Markov  random  fields,  random  circles, 
random  cylinders,  surface  roughness 


1«.  ABSTRACT  fConlinue  on  reverse  if  necessary  and  identify  by  block  number i 

The  aim  is  to  describe  the  exact  shapes  of  objects  that  were  meant  to  be  circles  or 
cylinders.  The  shapes  are  modeled  as  random  fields  whose  parameter  spaces  are  the  intended 
shapes.  A  specific  random  field  on  a  true  circle  is  introduced  via  exponential  smoothing 
of  a  random  fnoise,^  on  the  circle,  with  stationary  and  independent  increments.  The  result 
is  a  stationary  stochastically  continuous  random  field.  When  the  noise  is  the  white  noise, 
the  result  is  the  unique  continuous  stationary  Gaussian  random  field  on  the  circle.  A 
similar  definition  for  cylinders  yields  infinite  dimentional  Ornstein-^Uhlenbeck  processes. 


gnc  FILE  COB 


■f'  DTIC 

£  L.  E  CT  E 


SEP  1  5  1986 


20.  OISTRIBUTION/AVAILABILITY  of  abstract 
UNCLASSIFIED/UNLIMITEO  SAME  AS  RPT.  □  OTIC  USERS  □ 


22c  OFFICE  SYMBOL 


A.FOSR-TR*  86-0  65  7 


RANDOM  CIRCLES  AND  FIELDS  ON  CIRCLES* 


Erhan  tpinlar 


Approved  for  public  release  ; 
distribution  unlimited. 


J.  G.  Wang 


c  *  p  FO"  "F  C  7F  *  "2  C”  ~'"TF'lr,F’r  F  Z  jFAF  Cd  (  J) 
:r~ICL  OF  !?."?!  J. ’A  ITT  A  L  TO  OTIC 
'b.is  technical  report  his  been  reviewed  eni  ig 
eoproved  for  publ  ic  re  lease  IAA  A.-  R  190-12. 

0 ’ st r ibut i on  i s  uni imit ed. 

•'TTHE'A  J.  KE^T? 

'hief ,  Techrio  il  Inf  oraiation  Division 


Princeton  University 
Department  of  Civil  Engineering 
Princeton,  New  Jersey  08544 


Accession  For 

NTIS  GKA&I 
DTIC  TAB 
Unannounced 
Just  if lcatioo. 


Engineering-Economic  Series 
Report  EES-86-7 

May,  1986 


By - 

Distribution/ 

Availability  Codes 
(Avail  and/or 
Dist  I  Special 


r—» — ; — i _ i _ 

♦Research  supported  by  AFOSR  Grant  No.  82-0189  to  Northwestern  University. 


QUALITY 

INSPECTLD 


«6  9  15  069 


RANDOM  CIRCLES  AND  FIELDS  ON  CIRCLES* 


by 

E.  CINLAR  and  J.  G.  WANG 

1.  INTRODUCTION 

The  overall  aim  is  to  describe  the  exact  shapes  of  objects  that  were 

meant  to  be  circles,  cylinders,  and  so  on.  The  motivation  came  from  the 

need  to  model  small  deviations,  from  the  intended  surfaces,  caused  by  random 

effects.  For  instance,  in  high  precision  gas  lubricated  bearings,  the 
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clearances  are  of  the  order  of  10  inches  and  small  discrepancies  of  the  order 
of  10  6  and  10  7  must  be  taken  into  account  for  computing  lubricant  pressures 
and  load  capacities.  Aside  from  such  practical  considerations,  we  are 
motivated  by  a  desire  to  make  precise  the  notion  of  a  random  circle,  whose 
realizations  would  be  called  simply  circles  by  the  proverbial  man  in  the 
street. 

Analytically,  we  model  the  shapes  as  random  fields  whose  parameter  spaces 
are  the  intended  shapes.  Thus,  for  instance,  a  random  circle  is  a  random 
field  on  a  true  circle,  which  has  enough  stationarity  and  continuity  to 
deserve  the  term  circle. 

•Research  supported  by  AFOSR  Grant  No.  82-0189  to  Northwestern  University. 


-2- 


Preliminaries 

Throughout,  C  will  denote  a  circle  of  radius  one,  and  m  will  be  the 
length  measure  on  it,  that  is,  m(A)  is  the  length  of  A  for  each  arc  A  C  C. 

For  p,  q  €  C,  we  write  (p,q)  for  the  open  arc  going  from  p  to  q  counterclockwise 

and  we  write  m(p,q)  for  the  corresponding  arc  length.  Note  that, 

(p,p)  =  C\{p}  and  m(p,p)  =  2tt  .  The  notations  [p,q]  ,  (p,q]  ,  [p,q]  ,  etc.  are 
self-explanatory  arcs.  We  take  a  point  of  C  and  distinguish  it  by  labeling 

it  0.  We  associate  the  point  p  6  C  with  the  arc  length  m(0,p) ,  thus 

identifying  C  with  [0,2tt],  points  0  and  2tt  being  the  same.  For  p,q  6  C,  we 

define  p+q  as  the  point  whose  numerical  value  is  m(0,p+q)  =  m(0,p)  +  m(0,q) 

modulo  2tt .  Finally,  we  write  -p  for  the  point  satisfying  (-p)  +  p  =  0. 

Any  function  f  on  C  can  be  regarded  as  a  function  on  I?  =  (-00  ,  00)  by 
identifying  C  with  [0,2tt]  and  then  extending  f  onto  Rby  periodicity,  i.e., 
by  setting  f(2Trn  +  x)  =  f  (x)  for  all  O  <  x  <_  2tt  and  all  integers  n.  Continuity, 
right-continuity,  etc.  for  a  function  (and  therefore  a  stochastic  process)  on 
C  is  to  be  understood  in  the  usual  sense  for  the  periodic  extension  of  the 
function  onto  F. 

Finally,  throughout,  (fi,  f,  P)  will  be  the  complete  probability  space 
that  all  probabilistic  terms  refer  to. 

Random  fields  and  circles 

A  random  field  on  C  is  a  real-valued  stochastic  process  having  C  as  its 
parameter  space.  Each  random  field  X  on  C  can  be  thought  as  a  stochastic 
process  on  1  that  is  periodic  with  period  2tt  (that  is,  X2TTn+t  =  Xt 
for  all  t  and  all  integers  n) .  With  this  identification,  continuity,  right- 
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continuity,  stochastic  continuity,  stationarity ,  and  reversibility  carry 
over  to  random  fields  on  C.  We  repeat  some  of  these  once  more  in  the 
following. 

(1.1)  DEFINITION.  Let  X  =  (X  )  _  be  a  random  field  on  C.  Then, 

p  pec 

a)  X  is  said  to  be  stationary  if  the  probability  law  of  (X  )  c 
is  the  same  for  all  q  e  C; 

b)  X  is  said  to  be  reversible  if  the  probability  law  of  (X  )  c 
is  the  same  as  that  of  X; 

c)  X  is  said  to  have  the  Markov  property  for  the  arc  [p,q]  if  the 

collections  (X  )  .  and  (X  )  .  ,  are  conditionally  independent  given 

r  re lp,q]  r  re lq,pl 

Xp  and  x^;  X  is  said  to  be  Markov  if  it  is  Markov  for  every  arc  [p,q]  C  C. 

(1.2)  DEFINITION.  A  random  circle  is  a  positive,  stationary,  and  stochastically 
continuous  random  field  on  C. 

Each  realization  of  a  random  field  on  C  is  a  curve  on  the  surface  of  a 
cylinder  with  base  C.  However,  in  the  case  of  random  circles  at  least,  the 
point  of  view  we  adopt  is  that  shown  in  the  figure  below,  namely,  each  realiza¬ 
tion  of  a  random  circle  is  the  boundary  of  a  star-shaped  set. 

The  stationarity  refers  to  the  invariance  of  the  probability  law  under 
the  rotations  of  C,  or  equivalently,  under  shifts  of  the  "origin"  0  on  C. 

The  reversibility  refers  to  the  sameness  of  the  probability  law  whether  the 
figure  is  viewed  from  the  front  or  the  back  of  the  page.  In  the  absence  of 
"past"  and  "future"  on  C,  the  Markov  properly  (1.1c)  is  the  most  appropriate. 


and  it  is  the  most  common  way  that  property  is  defined  for  random  fields. 
Namely,  for  each  arc,  the  values  in  the  interior  and  exterior  of  the  arc  are 
conditionally  independent  given  the  values  at  the  two  boundary  points. 

For  a  stochastic  process  defined  on  CHAY  (1972)  used  the  phrase 
"quasi-Markov  on  [0,T]"  to  mean  that  for  each  sub- interval  [a,b]  of  [0,T] , 
the  values  of  the  process  in  the  interior  of  [a,b]  and  the  exterior  [0,T]\[a,b] 
of  [a,b]  are  conditionally  independent  given  X^X^.  It  is  easy  to  show  that 
the  usual  Markov  property  (of  the  past-future  type)  implies  the  quasi-Markov 
property  on  all  of  E.  Also,  it  is  clear  that  every  random  field  that  is 
Markovian  on  C  can  be  extended  to  a  process  on  S  that  is  quasi-Markov  on 
[0, 2ir]  . 


corresponding  radius  of  the  random  circle 
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Examples  and  comments 

Let  W  be  a  Wiener  process  on  B+,  let  a  >  0  be  a  fixed  number,  and 

define 

U.3I  *t  '  wt+a  '  V  Xp-Yp,  pec. 

The  process  Y  on  J?+  is  continuous,  stationary,  Gaussian,  and  quasi-Markov 
(the  latter  property  was  noted  first  by  SLEPIAN  (1961)).  Of  course,  X  is 
Gaussian  and  Markovian  on  C.  But,  X  is  not  continuous  on  C,  having  a  jump 
at  O  with  probability  one,  and  of  course,  not  stationary  on  C. 

Stationarity  on  C  is  easy  to  achieve.  It  amounts  to  picking  the 
distinguished  point  O  at  random  on  C.  Let  Y  be  an  arbitrary  random  field  on 
C,  let  U  be  independent  of  Y  and  have  the  uniform  distribution  on  C,  and 
set 


(1.4)  Xp  =  Yp_y ,  pec  . 

A  quick  computation  with  characteristic  functions  shows  that,  then,  X  is 
stationary  on  C.  In  particular,  if  Y  is  positive  and  continuous,  then  X  is 
a  continuous  random  circle.  More  particularly,  let  Y  describe  the  largest 
square  contained  in  C  and  having  0  as  a  corner;  then,  X  is  a  continuous 
random  circle  according  to  Definition  (1.2),  even  though  every  realization 
of  X  is  a  perfect  square.  This  is  somewhat  unpleasant  but  seems  inescapable 
Our  mental  picture  corresponds  to  (1.2)  when  the  variance  is  small  coitpared 


with  the  mean  radius. 


Let  X  be  a  Gaussian  random  field  on  C.  In  order  for  X  to  be  stationary 


on  C,  we  must  have 

(1*5)  Cov (X  , X  )  =  g°m(p,q)  =  gcrm(q,p) 

P  q 

for  some  covariance  function  g.  This  implies,  in  particular,  that  every 
stationary  Gaussian  random  field  on  C  is  reversible.  For  non-Gaussian  fields, 
stationarity  implies  (1.5)  but  (1.5)  does  not  imply  reversibility  (since 
covariances  no  longer  determine  the  probability  law) .  For  instance,  take  Y 
in  (1.4)  to  be  increasing  in  the  counterclockwise  direction  on  (0,2n)  and 
note  that  X  is  not  reversible. 

The  form  of  stationary  Gaussian  random  fields  on  C  is  well-known. 

LEVY  (1951)  has  shown  that  every  such  field  X  with  mean  zero  can  be  expressed 
as 


(1.6) 


X 

p 


00 


l 


n=0 


c  [U  cos  np  +  V 
n  n  n 


sin  np] , 


pec. 


where  the  U  and  V  are  all  independent  Gaussian  random  variables  with  mean  0 
n  n 

and  variance  1,  and  the  constants  c  are  square  summable. 

n 

Interestingly  enough,  each  term  of  the  series  (1.6)  has  the  Markov 
property  (1.1c)  for  most  arcs.  Consider 


(1.7)  Xp  =  U  cos  np  +  V  sin  np,  p^C, 

where  the  integer  n  >  1  is  fixed,  and  U  and  V  are  independent,  Gaussian,  with 
mean  0  and  variance  1.  Clearly,  X  is  continuous,  stationary,  Gaussian. 


It  has  the  Markov  property  for  every  arc  [p,q]  such  that  X  and  X  determine 

P  Q 

U  and  V,  which  is  whenever  sin  n(p-q)  /  0.  Otherwise,  if  np-nq  is  a  multiple 

of  tt,  the  Markov  property  ^or  [p,q]  fails.  For  instance,  if  n  =  1  and  p  =  TT /4 

and  q  =  5^/4,  then  cos  np  =  sin  np  and  cos  nq  =  sin  nq,  and  X^  and  X^  determine 

only  the  sum  U  +  V,  and  it  is  easy  to  see  that  the  Markov  property  fails. 

So,  this  X  is  not  Markovian  on  C. 

We  shall  show  that  the  process  X  given  by  (1.6)  with  c  =  l//2a 

o 

f~2  2~ 

c^  =  l//a  +  n  for  some  a  is  the  unique  Markovian  field  among  all  continuous, 
stationary,  Gaussian  ones,  up  to  multiplication  by  and  addition  of  constants. 

Integrators  over  C 

This  is  to  introduce  the  notation  and  terminology  we  shall  employ  with 
stochastic  integrals  with  respect  to,  essentially,  processes  with  stationary 
and  independent  increments. 

Throughout,  W  will  denote  the  white  noise  and  M  an  arbitrary  integrator 
with  stationary  and  independent  increments.  For  every  bounded  Borel  function 
f  on  C,  we  write 

(1.8)  Mf  =  /  f(p)  M (dp) 

C 

for  the  stochastic  integral  of  f  with  respect  to  M:  this  is  the  integral  of 
f  dZ  over  (0,2tt],  where  dZ^  =  M(dp)  and  Z  is  a  process  with  stationary  and 
independent  increments.  In  the  important  special  case  where  M  =  W,  the 


corresponding  Z  is  the  Wiener  process,  which  we  denote  by  (W  ) . 


(1.9)  REMARK.  Independence  of  increments  for  M  is  equivalent  to  saying 

that  Mf,  , ...,Mf  are  independent  whenever  the  f,  have  disjoint  supports. 
In  k 

Stationarity  for  M  is  equivalent  to  saying  that,  for  each  f, 

/  f (p)  M(dp) ,  J  f (p)  M (q  +  dp) 

C  C 

have  the  same  distribution  for  all  qEC.  We  express  this  by  the  phrase 
"the  probability  law  of  M (q  +  •)  is  free  of  q." 

The  structure  of  such  M  is  well-known:  for  every  bounded  Borel  function 
f  on  C , 


(1.10)  Mf  =  q.  mf  +  g  Wf  +  J  [N (dp,dz)  -  m(dp) n (dz) I  (z) ]  f (p)  z, 

CxK  (-1/-U 

where  a  and  g  are  constants,  m  is  the  arc  length  measure  and  W  is  the  white 

noise  as  mentioned  before,  and  N  is  a  Poisson  random  measure  on  Cx it  with  mean 

measure  m  x  n,  the  measure  n  being  a  Levy  measure  on  S  (that  is,  n-integral  of 
2 

z  — >  z  A  1  is  finite)  . 

Organization 

In  the  next  section,  we  show  that  the  random  field 

(1.11)  X  =  /  e~an,(p,q)  M (dq)  ,  p€C, 

P  C 

is  stationary,  right-continuous,  stochastically  continuous,  and  Markovian. 

In  the  special  case  M  =  W,  it  is  further  Gaussian,  continuous,  and  reversible. 
Conversely,  every  continuous  stationary  Gaussian  Markov  random  field  on  C 
has  the  form  c  +  bX  for  some  constants  c  and  b. 


If  M  is  positive  (which  makes  M  into  a  positive  compound  Poisson 
random  measure,  basically),  then  X  is  positive  and  is  a  Markovian  random 
circle.  In  the  case  M  =  W,  X  is  not  a  random  circle  because  it  does  take 
negative  values;  however,  for  small  enough  b  and  large  enough  c,  c  +  bX 
should  be  a  good  approximate  model  for  a  random  circle  with  mean  radius  c. 

In  general,  a  whole  class  of  Markovian  random  circles  Y  can  be  defined 
by  letting  Y  =  f(X)  for  continuous,  strictly  increasing,  strictly  positive  f 

In  section  3  we  concentrate  on  the  special  case  of  X  defined  by  (1.11) 
with  M  =  W,  the  white  noise.  The  process  X  satisfies  the  same  stochastic 
differential  equation  as  does  an  Ornstein-Uhlenbeck  process,  but  with  an 
unusual  boundary  condition,  namely,  that  the  initial  and  final  values  must 
be  the  same.  Therefore,  the  traditional  definitions  of  a  solution  (with 
the  initial  value  independent  of  the  Wiener  process )  do  not  apply.  We  do 
obtain  differential  equations  whose  solutions  in  the  traditional  sense 
yield  X.  Mere  interestingly,  we  give  a  decomposition 

\  =  f(f:!  xo  +  V  0  1  p  £  2  ' 

vnere  Y  is  a  Gaussian  process  independent  of  X^  and  is  Markovian  in  the 
ordinary  sense  (past-future  type) .  Finally,  in  the  same  section  we  give 
the  ourier  representat ion  for  X. 

Section  4  is  devoted  to  extending  X  onto  cylinders.  In  fact,  this  is 
not  difficult.  In  the  special  case  of  Gaussian  continuous  fields  on  the 
cylinder  C  >  1  we  obtain  an  Ornstein-Uhlenbeck  process  (Xt)  where  each  X^ 
is  a  random  field  p  - X  (p)  of  the  type  (1.11). 


The  main  difficulty  in  working  on  C  is  the  absence  of  past  and  future. 

We  have  been  able  to  avoid  the  issue  by  limiting  ourselves  to  integrands 
that  are  deterministic.  A  general  theory  of  stochastic  integration  on  circles 
is  yet  to  be  developed.  What  we  have  here  has  the  same  relation  to  that 
future  theory  as  does  the  Wiener  integrals  to  Ito  integrals. 


2.  STATIONARY  MARKOV  RANDOM  FIELDS  ON  C 

Throughout  this  section,  M  will  denote  an  integrator  on  C  with 
stationary  and  independent  increments .  The  particular  case  where  M  =  W, 
the  white  noise,  is  of  special  interest.  Recall  that  m  is  the  length 
measure  on  C.  Define 

(2.1)  X  —  f  e"am(p'q)  M(dq),  p€C, 

P  C 

where  a  ^  0  is  a  constant.  The  case  a  =  0  is  too  trivial  to  be  of  interest 
but  is  useful  in  rounding  out  the  theory.  The  following  is  the  main  result 
of  this  section. 

(2.2)  THEOREM 

a)  The  random  field  X  is  stationary,  right-continuous,  stochastically 
continuous,  and  Markovian  on  C. 

b)  In  the  special  case  where  M  =  W,  X  is  continuous,  Gaussian, 
and  reversible  in  addition  to  being  stationary  and  Markovian. 

c)  Conversely,  every  continuous,  stationary,  Gaussian,  and  Markovian 
random  field  Y  on  C  has  the  form  Y  =  b  +  cx,  where  b  and  c  are  constants  and 
X  is  as  in  (2.1)  with  M  =  W. 

If  it  were  positive,  X  would  be  a  Markovian  random  circle.  This  can  be 
achieved  by  taking  M  positive,  that  is,  by  letting  M  have  the  form 


(2.3) 


Mf  =  a 
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•  mf  +  j  N(dp,dz)  f (p) 2 
CxJi+ 

with  a  >_  0  and  the  Poisson  random  measure  N  on  having  the  mean  measure 
m  x  n  such  that  the  n- integral  of  z  A  i  is  finite.  We  state  this  next. 

(2.4)  COROLLARY.  Suppose  M  is  positive.  Then,  X  is  a  right-continuous 
Markovian  random  circle. 


The  discontinuities  of  X  coincide  with  the  atoms  of  M:  for  any 


p - >■  X^(oj)  ^as  a  jumP  size  -(1  -  e  27Ta)  z  at  q  if  and  only  if  the  measure 


N(co,*)  has  an  atom  at  the  point  (q,z).  By  taking  the  measure  n  vanishing 


-2ira. 


outside  an  interval  ic  K,  one  can  limit  the  jumps  of  X  to  sizes  -(1-e  )z 


with  zei.  Thus,  although  a  discontinuous  circle  is  repugnant,  discontinuities 
-6 


of  magnitude  10  are  not  likely  to  cause  comment  in  practice. 

However,  if  continuity  is  needed,  it  can  be  achieved  as  follows.  Let 
f:  » - >■  K+ be  a  continuous  strictly  increasing  function  and  define 


(2.5)  Y  =  f (X  ) ,  p€C  . 

P  P 

(2.6)  COROLLARY.  The  random  field  Y  is  a  right-continuous  Markovian  random 
circle.  In  particular,  if  M  =  W,  then  Y  is  also  continuous  and  reversible. 


Even  in  the  continuous  case,  by  choosing  f  appropriately,  one  can  limit 
the  values  of  Y  to  any  preselected  positive  open  interval  desired,  and  further 
one  can  make  the  distribution  of  Y  (necessarily  free  of  p  because  of 


stationarity)  any  continuous  strictly  increasing  distribution  function  one 

desires:  Choose  f  =  hog  where  g:  p  - >-  (0,1)  is  the  (cumulative)  Gaussian 

distribution  function  with  mean  0  and  variance  that  of  X  ,  and  where  h  is 

P 

the  functional  inverse  of  the  distribution  function  desired  for  Y  . 

P 

(2.7)  CONJECTURE.  Every  continuous  Markovian  random  circle  Y  has  the  form 
(2.5)  where  X  is  defined  by  (2.1)  with  M  =  W. 

If  true,  this  would  be  the  equivalent,  for  Markov  random  fields  on  the 
circle,  of  FELLER's  characterization  of  continuous  Markov  processes  on  P+. 

We  are  unable  to  prove  this,  but  we  suspect  strongly  that  the  conjecture 
is  very  near  the  truth. 

The  remainder  of  this  section  is  devoted  to  proofs,  basically. 

Stationarity 

This  is  easy  to  show  for  X.  The  following  does  it  for  a  larger  class. 

(2.8)  LEMMA.  Let  g  be  a  bounded  Borel  function  on  [ 0 , 2tt ]  and  define 


(2.9)  Z  =  J  g<?m(p,q)  M(dq)  ,  p6C  . 

P  n 


Then,  the  random  field  Z  is  stationary  on  C. 


PROOF.  Fix  sec,  and  note  that 


(2.10) 


2 


s+p 


J  g<->m(s+p,q)  M(dq) 
C 


=  f  g0m(s  +  p,  s  +  q)  M(s  +  dq)  =  J  gom(p,q)  M(s  +  dq) , 

C  C 

that  is,  the  random  field  Z  is  obtained  from  M(s  +  •)  by  the  same  rule  as 

s+* 

Z  is  from  M.  By  the  stationarity  of  M,  M(s  +  •)  has  the  same  law  as  M.  So, 

Z  has  the  same  law  as  Z. 
s+  • 


-ax 


The  random  field  X  has  the  form  (2.9)  with  g (x)  =  e  .  Another  special 
case  is  where  g(x)  =  x  and  M  =  W,  in  which  case 


(2.11) 


Z  =  Z0  +  2ttw  -  pW 


2v ' 


2Tr 

Z  0  =  J  s  dW  , 
U  0  S 


where  (W^)  is  a  Wiener  process,  in  other  words,  Z  is  a  tied  down  Brownian 
motion  with  an  initial  value  ZQ  that  depends  on  the  whole  of  W.  This  process 
was  introduced  by  LEVY  (1980,  vol .  5,  pp.  157  ff)  who  exploited  its  stationarity 
to  handle  the  Wiener  process  by  Fourier  transforms. 


Continuity 

This  is  easier  to  discuss  by  first  writing  the  integral  (2.1)  in  the 
more  traditional  form.  Identifying  C  with  (0,2tt),  separating  the  integral 
over  C  into  integrals  over  (0,p]  and  (p,2ir)  =  (0, 2ir]\(0,p]  ,  and  recalling 
that  m(p,q)  is  the  length  of  the  open  arc  (p,q)  going  from  p  to  q  counter¬ 


clockwise,  we  obtain,  for  p€(0,2^) 
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(2.12) 


(0,p3 


ap  -2TTa 
e  (e 


=  aP 


;  2TT-p+q) 

M(dq)  + 

J  e~a (q-p) 

<P,2tt] 

I  + 

/  - 

-  /  )  e'aq  M  (i 

(0,p) 

(0,2TT] 

(0,p] 

sap  (1  - 

-2na) 

e 

J  e  aq  M(dq) , 

(0,p] 


where,  as  in  (2.1), 


(2.13) 


X0  = 


e  aq  M(dq) 


(0, 2tt] 


The  last  integral  in  (2.12)  is  that  of  a  continuous  function  with 
respect  to  a  left-limited  right-continuous  process  over  (0,p] .  It  is  now 
obvious  that  X  is  left-limited  and  right-continuous  over  (0,2tt).  Taking 
limits  as  p  (  O,  we  see  that  X  is  right-continuous  at  p  =  0,  and  letting 
p  t  2tt,  we  see  that  X  has  a  left-limit  at  p  =  0  =  2Tr .  So,  X  is  left-limited 
and  right-continuous  on  C. 

It  follows  from  (2.12)  again  that,  for  every  outcome  weft, the  function 
P 


Xp(w)  has  a  jump  of  magnitude  -(1  -  e  2TTa)z  at  point  q  if  and  only 


if  the  measure  N(w,‘)  has  an  atom  (q,z) 
Thus,  p  - 


Xp(w)  is  continuous  if  and  only  if  N(w,*)  has  no  atoms. 


So,  p  — >  Xp(w)  is  continuous  for  almost  every  w  if  and  only  if  the  Levy 


measure  n  vanishes,  in  which  case  the  representation  (1.10)  for  M  becomes 
M  *  bm  +  cW  for  some  constants  b  and  c.  In  particular,  if  M  =  W,  X  is 


continuous . 


Finally,  since  the  mean  measure  of  N  is  m  x  n  and  m  is  diffuse,  for 
each  p€C,  the  probability  is  zero  that  N  has  an  atom  on  the  ray  {p}  x  F. 
Thus,  for  every  fixed  p€C,  X  is  continuous  at  p  almost  surely.  In  other 
words,  X  is  stochastically  continuous. 


Markov  property 


For  each  arc  (p,q]  C  C,  let  F  denote  the  a_algebra  generated  by 

Pg 

the  increments  M(A)  ,  A  C  (P/q)  »  and  let 


(2.14)  Z  =  /  e-am(P's)  M(ds)  . 

P<3  (p,ql 

Independence  of  the  increments  of  M  implies  that  F^  and  F^  are  independent. 

If  a  =  0  in  (2.1),  then  X  =  Xrt  for  all  p,  and  X  is  Markovian  on  C 

p  0 

trivially.  For  the  remainder  of  the  proof  we  take  a  ^  0. 

For  p  =  q,  the  Markov  property  for  [p,q]  is  trivially  true.  From  here 
on,  we  fix  p,qGC,  p  q. 

Writing  the  integral  (2.1)  as  the  sum  of  two  integrals,  one  over  (p,ql 
and  the  other  over  (q,p] ,  we  obtain 


(2.15) 


X 

P 


z 

pq 


+  e-am(p,q) 


X 

q 


e~am(q,p) 


Z  +  Z 

pq  qp 


Since  a  4  0,  Z  and  Z  are  determined  uniquely  by  X  and  X  through  these 

T  pq  qp  pq 

equations.  Hence,  the  a-algebra  generated  by  X^  and  X^  is  the  same  as  that 


generated  by  Z  and  Z 

pq  qp 


Therefore,  since  F  and  F  are  independent,  for  every  bounded 

pq  qp 


F  -  measurable  random  variable  F, 

pq 


E[F|X  ,  X  ,  F  ]  =  E[F|Z  ,  Z  ,  F  ]  =  E[F|Z  1 
1  p  q  qp  I  pq  qp  qp  1  pq 


and 


E  [F  |  X  ,X  ]  =  E  [F|  Z  ,  Z  ]  =  E  [ F I  Z  ]  , 

i  p  q  1  pq  qp  1  pq 


which  shows  that 


(2.16) 


E[F[X  ,  X  ,  F  ]  =  E[F|X  ,  X  ]  . 

1  p  q  qp  1  p  q 


In  other  words,  F  and  F  are  conditionally  independent  given  X  and  X  . 

pq  qp  p  ^ 


For  r  e  (p,q) ,  again  separating  the  integral  (2.1)  for  X^  into  integrals 


over  (p,q]  and  (q,p]  and  solving  (2.15)  for  Z  ,  we  obtain 

qp 


(2.17)  x  -  /  <Ta,"(r's>  „<dE>  ♦  z 


(p.qi 


qp 


/  M(dE)  +  _i_  (e-«"<r,q>  x  .  e-a»»(r,p)  , 

(P,q3  1-e  a  q  p 


It  is  evident  from  this  that 


(2.18) 


0(x  :  r  e  [p,q])C  0(X  ,  X  ,  F  ) 

X  F  4  r4 


and,  by  symmetry, 
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(2.19) 


0 (X  :  r  £  tp,q] )  C  0  (X  ,  X  ,  F  )  . 
r  p  q  qp 


In  view  of  (2.16),  namely,  because  the  o-algebras  on  the  right  sides  of 
(2.18)  and  (2.19)  are  conditionally  independent  given  X^  and  x^,  we  see 
that  the  desired  Markov  property  holds  for  the  arc  [p,q] . 


Completion  of  the  proof  of  Theorem  (2.2) 

The  preceding  sub-sections  have  shown  the  truth  of  (a) . 

In  the  special  case  where  M  =  W,  we  have  shown  that  X  is  continuous 
on  C.  Obviously,  in  addition,  X  is  Gaussian,  and  reversibility  is  immediate 
because  every  stationary  Gaussian  field  on  C  is  reversible.  So,  statement 
(b)  holds. 

To  show  (c) ,  let  Y  be  continuous,  Gaussian,  stationary,  and  Markovian 
on  C.  Subtracting  a  constant  and  dividing  by  another,  we  may  assume  that  Y 
has  mean  0  and  variance  1.  We  need  to  show  that  the  probability  law  of  Y 
is  that  of  bX  for  some  constant  b,  which  comes  to  showing  that  (since  both  X 
and  Y  are  Gaussian) , 


(2.20) 


EYY  =  b  EXX  ,  p  f  q  . 

p  q  P  q  T 


In  preparation  for  this  we  compute  the  covariance  of  X  and  X  from 

P  3 


(2.1)  with  M  *  W  and  a  ^  0.  We  get,  for  p  ^  q, 


(2.21) 


EX  X  =  E  f  f 
DO  J  J 


-am(p,r)  -am(q,s) 


W (dr)  W (ds) 


_  j  e~am (p, r) -am(q, r) 


m(dr) 


w 
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1_ 

2a 


(1 


-2TTa  -am(p,q)  -am(q,p) 
e  )  (e  +  e  )  , 


and 


(2.22) 


E(X  )2  =  (1  -  e'4^a) 

p  2a 


If  a  =  0,  we  have  X  =  W._  for  all  p,  and  hence  EX  X  =  2^.  Thus,  we  need 

p  2TT  p  q 

to  show  that,  for  some  a  >  0, 


(2.23)  EY  Y  -  (1  +  s-21-8)-1  (e-a"(p''3)  ♦  e-a"(q-Pl)  , 

p  q 


or  equivalently,  that  the  covariance  function  B,  defined  so  that  EY^Y^ 


B(m(p,q)),  has  the  form 


(2.24)  B (t)  =  be~at  +  (1  -  b)eat  ,  b  =  (1  +  e  2Tra)_1  . 


The  remainder  of  the  proof  follows  from  a  lemma  of  CHAY  (1972) .  Suppose 
that  Z  is  a  stationary  Gaussian  process  with  mean  0  and  variance  1  (stationarity 
over  H.)  .  Suppose  that  its  covariance  function  is  continuous.  Suppose  Z  is 
quasi-Markov  on  [0,T] .  Then,  the  covariance  function  f  of  Z  satisfies 
f " (t)  *  c  •  f (t)  for  t  £  T.  Depending  on  whether  c  =  0,  c  <  0,  or  c  >  0, 
there  are  three  possibilities  for  f (t)  : 


(2.25) 


1  -  ext,  cos  at, 


-at 

Ae  + 


(1  -  A)  e 


+at 
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where  a  >  0  in  the  first  two  cases  and  a  >  0  in  the  last  case  (and  there 
are  other  conditions  on  a,A,T). 

In  the  first  case,  the  corresponding  process  has  the  form  ^t+a  ~  W 
for  a  =  1/a.  This  process  is  not  continuous  on  C,  hence  1  -  at  cannot  be 
associated  to  Y. 

In  the  second  case,  cos  at,  the  corresponding  process  has  the  form 
U  cos  at  +  V  sin  at.  To  be  continuous  on  C,  a  >  0  must  be  an  integer. 

But,  then,  this  becomes  the  random  field  (1.7),  which  we  have  already  shown 
to  fail  the  Markov  property.  So,  the  covariance  function  of  Y  cannot  have 
the  form  cos  at. 

The  last  possible  form  is  in  fact  the  form  (2.24),  because  the  condition 
of  continuity  for  the  field  forces  the  values  at  t  =  0  and  t  =  2tt  to  be  the 
same,  which  condition  becomes 


(2.26) 


1  =  A  e 


(1  -  A) 


which  gives  A  =  b  and  a  =  a  as  in  (2.24) . 

This  completes  the  proof  of  Theorem  (2.2). 


Other  remarks 

Corollary  (2.4)  is  immediate  from  Theorem  (2.2a)  and  Definition  (1.2)  of 
random  circle.  The  remarks  on  discontinuities  of  X  following  Corollary 
(2.4)  were  proved  in  the  sub-section  above  entitled  "continuity . "  Corollary 
(2.6)  is  immediate  from  Theorem  (2.2)  and  the  following  facts.  First,  if  Z 
is  stationary,  then  so  is  f (Z)  for  any  field  Z  and  measurable  f.  Second, 


-2] 


similarly,  if  2  is  reversible,  so  is  f(Z).  If  Z  is  continuous  and  f  is 

continuous,  so  is  f  (Z) .  If  Z  is  Markovian  and  f  is  one-to-one,  then  f (Z) 

is  Markovian,  by  the  simple  fact  that,  then,  knowing  f (Z^)  is  the  same 

as  knowing  Z  and  vice-versa.  When  f  is  continuous  and  strictly  increasing, 
P 

f  is  obviously  one-to-one  from  B  onto  some  open  interval. 

Relationship  to  Ornstein-Uhlenbeck  process 

It  is  well-known  that,  when  is  the  parameter  space,  the  only 
stationary,  continuous,  Gaussian  Markov  processes  are  Ornstein-Uhlenbeck 
processes.  Since  the  random  field  X  has  the  analogous  properties  on  C  when 
M  =  W,  it  is  expected  that  X  be  the  analog  of  the  Ornstein-Uhlenbeck  process 
on  the  circle  C. 

Indeed,  the  representation  (2.12)  shows  this  clearly.  For  greater 
clarity,  we  re-write  (2.12)  as  a  stochastic  differential  equation  and  with 
M  =  W: 

(2.27)  dX  =  a  X  dp  -  b  dW  ,  0  <  p  <  2tt  , 

P  P  P 

_  2  'jjci 

where  b  =  1  -  e  .  (Of  course,  b  can  be  taken  to  be  arbitrary;  then, 
corresponding  solution  will  be  a  multiple  of  X.)  This  is  exactly  the 
stochastic  differential  equation  satisfied  by  Ornstein-Uhlenbeck  processes. 
But,  X  is  not  an  Ornstein-Uhlenbeck  process,  because  the  initial  random 
variable  XQ  depends  on  M  =  W  through  (2.13).  Instead,  X  is  the  unique  solu¬ 
tion  of  (2.27)  satisfying  the  two-sided  boundary  condition 


1  V  V  *«*  v  V  * -P  V  V  ’  •  "  -  *>  *.»  „ 
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(2.28) 


X0  = 


2  n 


In  the  theory  of  stochastic  differential  equations,  it  has  been 
customary  (see  IKEDA  and  WATANABE  (1981)  for  instance)  tc.  make  the 
convention  that  the  initial  value  is  independent  of  the  driving  process 
W  or  M,  or  whatever.  On  the  circle,  the  initial  value  is  also  the  final 
value,  and  such  conventions  will  have  to  be  altered. 


Characteristic  functional  of  X 

We  end  this  section  with  a  formula:  for  any  finite  measure  X  on  C, 


(2.29)  E  exp  i  J  Xdp)  X  =  exp  J  m(dq)  g(J  X(dp)  e 
C  P  C  C 

where  g  is  the  exponent  function  corresponding  to  M,  that  is,  if  M  is 
represented  by  (1.10), 


(2.30) 


i(6)  =  iae  -  j  62  e2  +  /  n(dz)  [e1  -  1  -  i6z  I(  x  (z)]  . 


Taking  X  =  A  6  + 


u  P  r 


+  X  6  ,  where  p.,  ...,  p  e  C  are  fixed  and 

n  p  0  n 

n 


X„,  ...»  X  G  D  are  constants  and  where  6  is  the  Dirac  measure  putting 
on  p 


its  unit  mass  at  the  point  p,  the  left-side  of  (2.29)  becomes  the 


characteristic  function  of  the  random  vector  (X  ,  ...,  X  ).  In  particular, 


taking  X  =  X  6  with  p  fixed,  we  get 
o  p 


E  exp(i  XQ  X  ) 


=  exp  j  m(dq)  g(XQ  e 


-am (p,q) 


C 

2TT 


=  exp  /  m (dq)  g(X 
n  u 


e-aq) 


(2.31) 


Two  special  cases  of  M  deserve  special  mention:  if  M  =  W,  the  white  noise 
2 

we  have  g ( 0 )  =  -}6  ;  and  if  M  is  positive  like  in  (2.3),  we  have 

00 

(2.32)  g(8)  =  ia9  +  /  n(dz)  (e1®2  -  1)  . 

0 

To  show  (2.29),  first,  we  observe  that 

/  X(dp)  X  =  /  M(dq)  /  X(dp)  ealtl(P,c3)  =  J  M(dq)  f  (q)  =  Mf 
with  an  obvious  definition  for  f.  Thus,  it  is  sufficient  to  prove  that 


(2.33) 


E  eiMf  =  em(<30f) 


for  every  bounded  Borel  function  f  on  C.  Various  continuity  arguments 
reduce  this  to  showing  that  (2.33)  holds  for  f  having  the  form 


f  =  6  1  +  ...  +  6  1  where  A.,  ...,  A  are  disjoint  arcs.  But,  then 

1  A,  n  A  1  n 

1  n 


iMf 


n 


E  e  =  A  E(exp  iS^  M (A^) ) 
k=l 


n 

=  H  exp(m(A^)  gtQ^))  =  exp  J  m(dq)  g ( f (q) ) 


as  required. 


3.  STANDARD  GAUSSIAN  FIELD 


This  section  is  devoted  to  the  random  field  X  defined  by  (2.1)  in  the 
special  case  where  M  =  W,  that  is, 

(3.1)  X  =  /  w (dq)  ,  p  €  C  . 

P  C 

We  exclude  the  trivial  case  a  =  0  and  assume  that  a  >  0.  The  random  field  X  is 
stationary,  reversible,  continuous,  Markovian,  and  Gaussian  with  mean  0  and 
covariance  (see  (2.21)) 


(3.2) 


B(p.q)  =  E  X^  X„  =  —(1  -  e 


p  q 


■2va^  (e-am(p,q) 


e-am(q,p) ) 


The  Markov  property  states  that,  given  XQ  and  X^,  the  processes 
(X^_:  r  €  [0,p])  and  {X^  8  [p,2tt])  are  conditionally  independent  given  X^.  Thus, 
given  X^,  the  conditional  law  of  X  is  the  law  of  an  ordinary  (past-future  type) 
Markov  process.  Our  first  aim  is  to  identify  that  ordinary  Markov  process. 

We  had  mentioned,  in  the  last  paragraphs  of  Section  2,  that  X  is  like 
an  Ornstein-Uhlenbeck  process,  but  it  is  not  one  because  the  stochastic  dif¬ 
ferential  equation  to  be  solved  is  driven  by  a  Wiener  process  W  that  depends  on 
the  initial  value  X^.  Our  second  aim  is  to  obtain  the  equations  for  X  whose 
driving  term  is  independent  of  X^. 

Finally,  we  shall  obtain  the  Fourier  representation  for  X,  which  turns  out 
to  be  almost  as  slow  as  that  of  the  Wiener  process. 


Before  listing  the  results,  we  introduce 

2tt  ci 

<3-3)  b  =  1  -  e  ;  =  a  m(p,q) ,  p,q  6  C  , 

in  order  to  lighten  the  notation  somewhat.  We  shall  be  using  hyperbolic 
functions  extensively. 


(3.4)  THEOREM.  Let  [p,q]  be  a  fixed  closed  arc  on  C  with  0  <_  p  <  q  <  2t . 

/N 

There  exists  a  Wiener  process  W  on  [p,q]  independent  of  X^  and  X^  such  that 


sinh  m 


(3.5) 


X  = 


sinh  m 


X  + 


X  + 


r  sinh  m  p  sinh  m  q 
pq  pq 


»  / 


sinh  m 


,  ,  sinh  m 

tP.r]  sq 


W (ds) 


for  all  r  €  [p,q] . 


(3.6)  COROLLARY.  Let  Y^  be  the  last  term  on  the  right  side  of  (3.5). 
Considered  as  a  process  on  the  time  interval  [p,q] ,  Y  is  a  Markov  process 
(in  the  ordinary  sense),  is  independent  of  X  and  X  ,  and  has  Y  =  Y  =  0. 

p  q  p  q 

Moreover,  it  satisfies  the  stochastic  differential  equation 


(3.7) 


dY  =  -a  (ctnh  m  ) 
r  rq 


Yr  dr  +  bW (dr) , 


p  <  r  <  q 


(3.8)  COROLLARY.  Given  X  and  X  ,  the  conditional  law  of  (X  )  .  , 

P  q  r  r€[p,q] 

law  of  an  ordinary  Markov  process  on  [p,q] .  The  process  satisfies  the 


is  the 


stochastic  differential  equation 


-26- 


(3.9)  dX  =  (aX  csch  m  -  aX  ctnh  m  )dr  +  bW(dr),  p  <  r  <  q  , 
r  q  rq  r  rq 


where  W  is  a  Wiener  process  independent  of  the  pair  (X^.X^) . 


(3.10)  PROPOSITION.  There  exists  a  Wiener  process  W  independent  of  XQ  such 
that 


,  .  sinh  m  _ 

»•“>  *-  -  **»£+**  k0  .  b  J  _  S.a,» , 


cosh  ita  0 


p  6  C 


[0,p]  s,2tt 


Thus,  given  XQ,  the  conditional  law  of  X  is  the  law  of  an  ordinary  Markov 
process  on  the  time  interval  [0,2iT]  . 


(3.12)  remark.  Still  other  facts  about  X  can  be  obtained  from  Corollaries 
(3.6)  and  (3.8)  by  replacing  p  and  q  there  with  0  and  2 tt  respectively. 

(3.13)  THEOREM.  The  random  field  X  has  the  series  representation 


(3.14)  X  =  — Un  +  — 

P  a/2:TT  JtT  n=l  V"  2  2 

/a  +n 


00 

_  1 

}  -  [U^  cos  np  +  V_  sin  np]  , 

_ .  r  ^  n  n 


where  Uq,U^ , . . . ,V^, , . . .  are  independent  Gaussian  random  variables  with 

mean  0  and  variance  1.  In  fact  Un  =  W„  //2tT,  and  for  n  4  0/ 

0  t  it  ' 


2tt 


(3.15)  U_  =  - - -  /  [a  cos  nq  +  n  sin  nq]  W(dq)  , 

/7r(a  +n  ) 


2tt 


V  = 
n 


- —  /  [a  sin  nq  -  n  cos  nq]  W(dq) 


7”  2,  0 


/rr(a  +n  ) 


The  remainder  of  this  section  is  devoted  to  proofs.  Corollaries  (3.6) 


and  (3.8)  are  easy  consequences  of  Theorem  (3.4).  Proposition  (3.10)  is 
obtained  from  Theorem  (3.4)  and  Corollary  (3.8)  by  replacing  p  with  0  and 
q  with  2tt.  Thus,  the  only  things  to  be  proved  are  (3.4),  the  Markov  representa 
tion  theorem,  and  (3.13),  the  Fourier  series  representation. 


Proof  of  Theorem  (3.4) 


For  r  G  [p, q] ,  define 


sinh  m 


(3.17) 


y  =  x  - 


-IR 


sinh  m 


X  - 


_EI 


r  r  sinh  m  p  sinh  m  q 

pq  pq 


(3.18) 


w  =  r  Y  +  f  /  Y  ctnh  m  du 
r  b  r  b  [p, r]  U  u* 


It  is  clear  that  Y  is  a  continuous  Gaussian  process  on  [p,q]  with  mean 

A 

0  and  Y  =  Y  =0.  So,  W  is  again  a  continuous  Gaussian  process  with  mean  0, 
P  ^3 

A 

and  =  0.  In  differential  form,  (3.18)  gives 


d  (—  -  Y  )  =  - -  (dY  + 

sinh  m  r  sinh  m 

ra  rq 


cosh  m 


.  ,  a  Y  dr) 

r  sinh  m  r 

rq 


— r-r -  dW 

sinh  m  r 
rq 


Hence,  since  Y  =  0, 
P 


Y  =  sinh  m  f  ■—  -  dW 

r  rq  .  1  ,  sinh  m  s 


(p.r) 


sq 


Putting  this  into  (3.17)  and  rearranging  the  terms  yield  (3.5). 
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To  complete  the  proof,  we  need  to  show  that  W  is  a  Wiener  process  and 
is  independent  of  X^  and  X^.  Since  all  the  processes  are  Gaussian,  these 
matters  are  reduced  to  computing  various  covariances. 

First  we  express  (3.2)  in  different  notation: 


(3.19) 


E  X  X  =  —  e  ra  cosh  (Tia  -  m  )  ,  r,s  6  C 
r  s  a  rs 


Using  this  with  (3.17)  and  some  elementary  identities  yields 


(3.20) 


(3.21) 


EYX  =  E  Y  X  =  0, 
r  p  r  q 


r  6  [p,q]  , 


E  Y  Y  =  - — - 

r  s  a  sinh  m 


sinh  m  sinh  m  ,  p<r<s<q 
pr  sq  —  — 


Since  X  and  Y  are  Gaussian,  (3.20)  implies  the  independence  of  the  process 
Y  from  the  pair  (X^,X^) ,  and  in  view  of  (3.18),  the  process  W  is  independent 
of  X  and  X  . 

P  q 

A  A 

Since  =  0  and  since  W  is  a  continuous  Gaussian  process  with  mean  0, 
in  order  to  show  that  W  is  a  Wiener  process,  it  is  sufficient  to  show  that 


(3.22) 


E  W  W  =  r  -  p  , 
r  s  r 


p  £  r  <  s  <_  q  . 


This  computation  requires  some  delicacy,  for  which  reason  we  provide  a  few 


details. 


It  turns  out  to  be  useful  to  introduce  the  notations 


(3.23) 


F  =  J  sinh  m  ctnh  m  dv  , 
i  PV  vq 

p,u 


-29- 


(3.24) 


aG 

u 


J 

lp,u] 


a  cosh  m 

vq 


dv  =  sinh  m 

pq 


sinh  iti  . 
uq 


It  follows  from  (3.18)  that,  for  r  <  s, 

,  .r 
b  =  Vs  +  a  J  du  Vs  Ctnh  "uq 
P 

r  s 

+  a  (J  du  +  /  du)  YuYr  ctnh 

P  r 

r  u  s 

+  a  /  du  (/  d v  +  J  dv)  Y^  ctnh  muq  ctnh 
p  P  u 

from  which,  taking  expectations  with  the  aid  of  (3.21)  and  using  the  nota 

tions  of  (3.23)  and  (3.24),  we  obtain 

(3.25)  a  sinh  E 


=  sinh  m  sinh  m  +  a  sinh  m  /  dF 
pr  sq  sq  u 


r  s 

+  a  sinh  m  J  dF  +  a  sinh  m  J  dG 
rq  1  u  pr  u 


+  a2  J  dGu  /  dFv  +  a2  /  dFu  /  dGv 
P  P  P  u 


=  sinh  m  , 
Pr 


(sinh  m  +  a(Gg  -  Fy)  ] 


+  a  [sinh  m  +  sinh  +  aG^  +  aGg]  Ff 


Note  that,  directly  from  (3.23)  and  (3.24), 


2a 1  J  Gu  dFu  =  2aFr  sinh  mpg  -2a  /  sinh  mpu  cosh  mug  du 


=  2aF  sinh  m  -a  /  (sinh  m  +  sinh(m  -  m  )]du 
r  pq  i  PS  pu  uq 


2aF  sinh  m  -  a  •  (r  -  p)  sinh  m 
r  pq  PS 


-  4  cosh (m  -  m  )  +  4  cosh  m 
2  pr  rq  2  pq 


Putting  this  into  (3.25)  and  simplifying,  we  obtain  (3.22),  thus  completing 
the  proof  of  Theorem  (3.4)  . 


Proof  of  Theorem  (3.13) 

Considering  X  as  complex  valued,  Fourier  series  for  X  is 


(3.26) 


X  =  I  AneinP  '  P  6  c 

p  "  n 


where 


(3.27) 


A  -  i-  J  e'inP  X  dp 
r.  2tt  i  P 


It  will  be  convenient  to  introduce,  in  addition  to  b  =  1  -  e 


(3.28) 


c  =  — -  ;  C  =  1  '  11  t  0 


K-k! 


n  2rr 


/  W (dq)  e"am(P'^ 


1  /  W  (dq)  -K-  e'inq 

2tt  ^  a- in 


b  (• 

- 2 — 2~  J  W(dq) (a  +  in)  (cos  nq 

2tt  (a  +n  )  C 


i  sin  nq) 


Putting  Uq  =  and  defining  and  for  n  f  0,  for  n  positive  and 


negative,  by  (3.15)  and  (3.16),  we  see  that 


A  =  be  U  :  A  =  -q-  be  (U  -  iV  )  ,  n  ^  0  . 

0  00  n2nn  n  r 


Putting  these  into  (3.26)  yields 


(3.29) 


Xp  =  bc0U0  +  2  L  °n  (Un  “  iVn)  (coS  np  +  1  sin  np) 


■  bcnUn  +  I  I  c  [<U  +  U  -  iV  -  iV  )  cos  np 


n  n  -n  n  -n 


+  (iU  -  iU  +  V  -  V  )  sin  np] 
n  -n  n  -n 


bcnU  +  —  l  c  (2U  cos  np  +  2V  sin  np) 
uuii  .nn  n 

n=l 


since  U  =  U  and  V  =  -V  for  all  n  >  1.  This  is  exactly  (3.14),  and 
“n  n  -n  ri 

the  only  thing  remaining  is  to  show  the  claims  about  the  nature  of  U  and  V 


-  v*  r\*  i  ^ j  r j  r j  r j1  r r\  tv  w  wv  ru  » w 
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It  is  obvious  that  they  are  Gaussian  and  have  mean  0.  The  remaining 
assertions  about  their  variances  and  independence  from  each  other  are  proved 
by  showing  that 

2  2 

EU  =  E  V  =1,  EUU  =  E  V  V  =  0,  n  4  m  , 
n  n  n  m  n  m  1 

and  that 


EU  V  =  0,  n  >  0  ,  m  >  0  . 
n  m  — 


This  requires  only  elementary  calculus  once  we  note  that 


2tt 

E  J  f  (p)  W  (dp)  /  g  (q)  W(dq)  =  J  f  (p)  g  (p)  dp  . 
C  C  0 


(3.30)  REMARK.  One  could  use  the  fact  that  X  is  stationary  and  periodic 

and  EX  X  =  B(p  -  q)  where 
P  q 


b2  00  2  2 

B  (t)  =  - -  +  T  b  c  cos  nt  ,  t  €  3R  , 

_  2  n 

2ira  n=l 


and  then  use  the  spectral  decomposition  of  stationary  Gaussian  processes  to 
arrive  at  (3.14).  The  proof  above  i.>  longer,  but  it  identifies  the  coef¬ 


ficients  U  and  V  . 
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4.  FIELDS  ON  CYLINDERS 

Extending  the  results  of  the  preceding  sections  to  random  fields 
defined  on  cylinders  presents  no  new  difficulties.  We  limit  ourselves, 
therefore,  to  introducing  the  random  field  of  interest  and  to  stating 
several  interesting  properties. 

Let  C  be  the  circle  of  radius  one  as  before  and  consider  the  cylinder 
CxIR.  For  p  and  q  on  C,  we  write  m(p,q)  for  the  length  of  the  arc  (p,q) 
going  from  p  to  q  counterclockwise,  just  as  before.  In  addition,  we  write 
m(s,t)  for  the  length  of  the  interval  (s,t)  C  )R»  although  not  permissible  by 
the  strict  rules  for  notation  making. 

Let  M  be  an  integrator  on  C  x  H  with  stationary  and  independent 
increments:  for  every  bounded  Borel  function  f  with  compact  support  in  C  x  H, 

(4.1)  Mf  =  a  J  m(dp)  dt  f(p,t)  +6  J  W(dp,dt)  f(p,t) 

+  J  [N (dp, dt,dz)  -  m(dp) dtn (dz) I  .  (z)  ]  f(p,t)z  , 

Cx*XR  t-IfU 

where  a  and  g  are  constants,  N  is  a  Poisson  random  measure  on  Cxp.XR  with  mean 
measure  element  m(dp)  dt  n(dz),  where  n  is  a  Levy  measure,  and  W  is  the  white 
noise  on  CxR.  The  last  may  be  thought  as  the  integrator  corresponding  to 
Brownian  sheet,  that  is. 


W  =  J  W(dq,ds)  ,  0  <  p  <  2tt,  0  <  t  < 

P'  (0,p]x(0,t] 


H  \  > 


(4.2) 


00 
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is  a  Brownian  sheet  on  [0,2tt]  x  b+,  which  is  a  Gaussian  process  with  mean 
0  and 

(4.3)  E  W  W  =  (p  A  q)  (t  A  u)  . 

p,t  q,u 

Integration  with  respect  to  W  was  considered  by  WONG  and  ZAKAI  (1974) , 
CAIROLI  and  WALSH  (1975),  and  others.  The  last  integral  in  (4.1)  is  in  fact 
simpler  and  should  be  understood  in  analogy  with  the  simple  time  case.  The 
integral  over  Cx  :R  x  (R\(-l ,  1) )  =  e  is  that  of  f(p,t)z  with  respect  to  a 
random  measure,  and  hence,  is  an  ordinary  integral.  The  integral  on  the 
complement  of  E  is  a  stochastic  integral;  it  is  the  limit  in  probability 
(indeed  almost  sure  limit)  of 

f  [N (dp,dt,dz)  -  m(dp)  dt  n(dz)]  f(p,t)z  , 

JEe 

as  C  i  0,  where  E£  =  C  *  K  * ( (-1 , 1)\(-E ,e) )  . 

The  random  field  of  interest  is 

(4.4)  X  (p,  t)  =  /  e  bm(s,t)  M(dq>cjsjf  p  €  C,  ten, 

Cx  (-00,  t] 

where  a  ^  0  and  b  >  0  are  fixed  constants.  We  view  it  as  a  random  field  on 
the  cylinder  Cx  k.  Alternately,  we  may  regard  it  as  the  evolution  in  time 
of  a  random  field  on  the  circle  C. 
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The  message  of  the  following  theorem  is  that  the  sections  of  X  are  all 
stationary  and  Markov.  Recall  that  a  time-homogeneous  Markov  process,  with 
time-set  E+  and  taking  values  in  a  Lusin  space,  is  called  a  Hunt  process  if 
it  is  right-continuous,  left-limited,  strong  Markov,  and  quasi-left-continuous, 
the  last  two  properties  being  with  respect  to  a  filtration  that  is  right- 
continuous  and  augmented  properly.  See  BLUMENTHAL  and  GETOOR  (1968)  for  the 
definitions.  We  shall  use  the  term  here  for  processes  with  time-set  E,  with¬ 
out  specifying  the  filtration,  and  somewhat  fraudulently:  Extension  of  the 
concepts  onto  His  no  problem;  we  may  take  as  filtration  the  (F^)^  ,  where 

Ft  =  vu,  F°  being  the  a-algebra  generated  by  {Mf:  f  continuous  with 
compact  support  contained  in  C  x(-«>,t]}  and  H  being  all  the  null-sets  of  the 
completion  of  the  fraud  is  that  we  have  only  one  probability  measure  P, 
instead  of  a  whole  collection,  one  for  each  starting  state,  of  which  we  have 
none . 

(4.5)  THEOREM,  a)  For  each  t®E,  p  — >  X(p,t)  is  a  stationary,  right-continuous 
stochastically  continuous,  Markovian  random  field  on  C;  the  law  of  p  — >  X(p,t) 
is  free  of  t. 

b)  For  each  p  6  C,  t  — X(p,t)  is  a  stationary  Hunt  process. 

c)  For  any  integer  n  _>  1  and  points  p^,  —  ,p^  on  C,  the  process 

t  — >  (X (p^ , t) , . . . , X (p^ , t) )  is  a  stationary  Hunt  process  with  state  space  En. 

d)  Let  X ( • , t)  denote  the  mapping  p  — >  X(p,t).  Then,  the  process 

t  — >  X ( • , t)  is  a  stationary  Hunt  process  taking  values  in  the  space  E  of  all 
right-continuous,  and  left-limited  functions  from  C  into  E,  topologized  by 
uniform  convergence  on  C. 
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(4.6)  COROLLARY.  Suppose  that  M  =  W,  the  white  noise  on  Cx  B.  Then,  in 
addition  to  the  properties  above,  X  is  a  continuous  stationary  Gaussian  random 
field  on  C  xR  with  mean  0  and 

(4.7)  E  X  (p,  t)  X  (q,u)  =  B(p,q)  e  b|U  ,  p,q  G  C,  t,u  G  H  . 

where  B  is  as  given  by  (3.2).  Moreover, 

a)  for  each  p  G  C,  t  — ^  X(p,t)  is  an  Ornstein-Uhlenbeck  process; 

b)  for  any  integer  n  _>  1  and  points  p^,...,p  6  C,  the  process 

t  — >•  (X(p,  t),...,X(p  ,t))  is  an  pn  -  valued  Ornstein-Uhlenbeck  process, 

J- ,  n 

which  is  a  weak  solution  of 


(4.8) 


dX1 (t)  *  -bX1(t)  dt  + 


7  c..WD(dt) 

j-1  13 


where  X^(t)  =  X(p^,t),  and  w\...,wn  are  independent  Wiener  processes,  and 
the  matrix  c  satisfies  cc'  =  B,  that  is. 


(4.9) 


n 

kI1  Cik  Cjk  =  B(Pi'Pj>  '  1  <  i,j  <  n  ; 


c)  the  process  t — >  X(»,t)  is  an  Ornstein-Uhlenbeck  process  taking 
values  in  the  space  <E  of  all  continuous  functions  from  C  to  R,  topologized 
by  uniform  convergence  on  C. 


(4.10)  REMARK.  An  Ornstein-Uhlenbeck  process  with  values  in  lor  3*°  is  a 
process  that  satisfies  an  equation  like  (4.8),  see  IKEDA  and  WATANABE  (1981) 
for  instance.  In  the  case  of  processes  taking  values  in  infinite  dimensional 


spaces  like  C  in  (4.6c),  see  ITO  (1982)  and  MEYER  (1982)  for  definitions. 


Our  case  is  much  clearer,  because  of  the  explicit  formula  (4.4). 

(4.11)  THEOREM.  The  random  field  X  enjoys  the  following  Markov  property. 

Let  A  Q  C  be  an  arc,  let  B  Q  P  be  an  interval,  and  consider  the  rectangle 

A x  B  on  the  surface  of  the  cylinder  Cx  J*.  Then,  {X(p,t):  p  6  A  and  t  6  B} 
and  { X (p , t ) :  p  0  A  or  t  g  B}  are  conditionally  independent  given  {x(p,t) : 
(p,t)6  3  ;A  xB)},  3  (A  x B)  being  the  boundary  of  Ax  B. 

(4.12)  REMARK.  The  arc  A  can  be  taken  to  be  C,  in  which  case  Ax  B  is  a 
cylinder,  and  the  Markov  property  above  coincides  with  the  statement  (4.5d). 
The  degenerate  arc  A  =  {p}  is  allowed.  Similarly,  B  does  not  have  to  be 

an  interval. 


(4.13)  REMARK.  By  limiting  the  parameter  t  to  an  interval  D  ^  I?,  we  may 
regard  X  as  a  field  on  the  cylinder  Cx  D.  For  this  purpose,  D  =  =  (0,<») 

and  D  =  [0,d] ,  d  >  0,  are  of  special  interest.  All  the  results  listed  above 
continue  to  hold  for  such  an  X.  In  considering  such  cylinders  CxD,  it  may 
appear  more  natural  to  define  X  via  integrators  on  CxD.  Assuming  D  =  [0,d] 
or  D  =  (0,o°)  ,  this  can  be  done  by  setting 


(4.14) 


X (p,t)  =  e  bt  X (p,0)  +  J 

Cx (0, t] 


a-.«(P,,)-b.(.,t)  M(aq-ds), 


where 


(4.15) 


X  (p,0)  =  /  e-am(P'<3)  M  (dq)  , 


Mq  being  an  integrator  on  C  independent  of  the  integrator  M  on  C><D,  Of 
course, 

0  bs 

(4.16)  Mq  (dq)  =J  e  M(dq,ds)  . 

—  00 

The  formula  (4.14)  can  be  used  to  define  an  "Ornstein-Uhlenbeck" 

process  t  - »X(»,t)  on  IR+ and  taking  values  in  E  of  (4.5d)  with  initial 

state  X  (*,0)  arbitrary  and  independent  of  M  on  C*B+.  In  an  obvious  sense 
the  probability  law  of  X(*,t)  approaches,  as  t  — >  °°,  to  that  of  the  right- 
side  of  (4.15)  with  Mq  having  the  same  law  as  the  right-side  of  (4.16). 

Many  of  the  results  above  rest  on,  or  motivated  by,  the  following 
observation.  Its  truth  is  immediate  from  Remark  (1.9)  and  a  similar  state¬ 
ment  for  integrators  on  JR. 

(4.17)  LEMMA.  Define  M(gxh)  by  (4.1)  with  f(p,t)  =  g  (p)  h(t),  where  g  is 
a  bounded  Borel  function  on  C  and  h  is  bounded  Borel  function  with  compact 
support  in  JR.  Then, 

a)  for  fixed  g,  h  — >  M(g*h)  is  an  integrator  on  JR  with  stationary 
and  independent  increments; 

b)  for  fixed  h,  g  — >  M(gxh)  is  an  integrator  on  C  with  stationary 
and  independent  increments. 

In  a  similar  vein,  and  more  interesting,  is  the  following  "stochastic 
kernel."  For  h  bounded,  Borel,  with  compact  support  in  JR,  define 


-39- 


(4.18)  K(p,h)  =  /  e"am(p'q)  •  h(t)  M(dq,dt)  . 

Cx  £ 

Recall  that  £  stands  for  the  space  of  all  right-continuous  left-limited 
functions  from  C  into  £  (see  (4.5d)).  The  following  is  immediate  from  the 
preceding  lemma . 

(4.19)  PROPOSITION.  For  each  fixed  h,  p  — >  K(p,h)  is  a  random  element  of 

I.  For  each  fixed  p,  h  >  K(p,h)  is  an  integrator  on  £  with  stationary 

and  independent  increments. 

Proof  of  Theorem  (4.5) 

Stationarity  of  the  increments  of  M  means  that  the  probability  law  of 

M (p  +  dq,  t  +  du)  is  free  of  (p,t) .  This  implies  that  the  law  of  X  is  free 

of  the  choice  of  origin  and  of  rotations  of  the  coordinate  system.  This 

shows  the  stationarity  in  all  the  statements  of  (4.5). 

For  fixed  t  6  B,  let  h(s)  =  e  for  s  <  t  and  h(s)  =0  for  s  >  t. 

Then,  by  Lemma  (4.17b) ,  g  — >  M(g  xh)  is  an  integrator  on  C  with  stationary 

and  independent  increments,  and  p  — >  X(p,t)  is  defined  from  it  according  to 

the  formula  (2.1).  Thus,  (4.5a)  is  a  re-statement  of  Theorem  (2.2a). 

For  fixed  p  6  C,  let  g(q)  =  e  am^P»^ .  By  Lemma  (4.17a),  h  — >  m  (h)  = 

P 

M(gxh)  is  an  integrator  on  £  with  stationary  and  independent  increments,  and 
we  have 

(4.20)  X(p,t)  =  /  e-fc>m(s,t)  M  (ds)  _ 


That  t 


X(p,t)  is  left-limited  and  right  continuous  is  immediate  from 


this.  We  have 


(4.21)  X (p, t+u)  =  e'bu  x(p,t)  +  /  e“bm(s'u)  Mp(t+ds)  . 

0 


Since  M^(t+«)  is  independent  of  the  past  and  has  the  same  law  as  M^(*), 
(4.21)  implies  that  t  — — >  X(p,t)  is  a  time -homogeneous  Markov  process. 

Strong  Markov  property  is  similar;  one  can  replace  t  by  a  stopping  time  T 
in  (4.21)  and  use  the  known  "strong  Markov  property"  for  M^.  Finally,  quasi¬ 
left  continuity  comes  from  the  facts  that  every  jump  of  t  — >  X(p,t)  coincides 

with  an  atom  (t,z)  of  the  Poisson  random  measure  involved  with  M  and  there- 

o 

fore,  the  jump  times  are  exhausted  by  a  sequence  of  stopping  times  each  one 
of  which  is  totally  unpredictable.  This  completes  the  proof  of  (4.5b). 

The  proofs  of  (4.5c,d)  are  entirely  similar;  To  show  (4.5d),  replace 
Mq  in  the  preceding  paragraph  with  h  — s.  K(*,h),  the  latter  being  defined  by 
(4.18)  and  having  the  properties  mentioned  in  Proposition  (4.19). 


Proof  of  Corollary  (4.6) 

When  M  =  W,  the  white  noise  on  cx  E,  it  is  obvious  that  X  becomes 
Gaussian,  stationary,  reversible,  continuous.  It  is  easy  to  compute  the  co- 
variance  function  and  show  that  (4.7)  is  indeed  true. 

The  remaining  statements  are  all  similar.  We  prove  (b)  only.  First, 
putting  M  *  W  and  taking  differentials  in  (4.20)  we  have 


(4.22) 


dX (p, t)  =  -bX (p,t) dt  +  W  (dt) 


where 


(4.23) 


W  (h)  =  /  e 

P  Cxi? 


-am(p,q) 


h(t)  W(dq,dt)  . 


Let  p_,...,p  e  C  be  fixed.  For  h,,...,h  bounded,  Borel,  with  bounded 
in  in 


support ,  we  have 


E  W  (h.)  W  (h.)  = 

p .  1  p  "i 

l 


*  /  exp[-am(p. ,q)  -  am(p.,q))  h.(t)  h.(t)  m(dq)dt 


C  x  B 


=  B(Pi,p  )  /  hi(t)  h  (t)  dt 


-E«j  cikwk<hi>>  <J.  Sk’ik(hi))  ■ 


k=l 


jk 


with(c. .)  and  )  as  described  in  (4.6b).  Thus,  the  probability  law  of  the 

k 


n-dimensional  process  (VJ  ,  —  ,W  )  is  that  of  (  T  c.,  W  ).  ,  .  Putting 

P1  pn  k=l  lk  1=1* ...»n 

this  fact  together  with  (4.22)  we  obtain  (4.8)  and  the  proof  of  (4.6b)  is 


complete . 


Proof  of  the  Markov  property  (4.11) 


It  is  sufficient  and  convenient  to  give  the  proof  for  A  =  [0,p]  ^  C 
and  B  *  |0,d]C»,  d  >  0.  We  let 

(4.24)  L  =  C  x  (-oo,0]  ,  G  =  [0,p]  x  [0,d]  , 

H  =  [p,0]  x  [0,d]  ,  K  =  Cx  [d,co)  . 
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r*4  .’I 


Ji  Ak  Jl  .•!  . 


For  any  collection  {2{t):  t  6  C x  B}  and  any  set  SCCx  B,  we  write  Z  for  the 

b 

o-algebra  generated  by  {Z(t):  t  6  S} .  Similarly,  we  write  M  for  the  a-algebra 

b 

generated  by  {Mf:  f  =  ly,UCS}.  Finally,  if  E,  E,  X  are  sub-o-algebras  of  F, 
we  write  (read  E  splits  E  and  X) 


E  ]  E  [  X 


to  mean  "E  and  X  are  conditionally  independent  given  E."  In  the  special  case 


when  E  and  X  are  independent,  we  write 


E  ]  [  X  . 


The  following  two  elementary  facts  will  be  used  a  number  of  times: 


(4.25)  E  v  B  J  E  [  X  =*  EvB]EvB[X=*  E  vB]EvB[XvB 


(4.26)  E  ]  E  [  X,  E  ]  [  E  v  E  V  X  E  ]  E  [  X  v  E 


We  start  by  letting 


(4.27) 


Y(p,t)  =  / 

Cx  (0,t] 


e-am(P,q)-bm(s,t)  M(dq;ds)  . 


Applying  the  arguments  of  the  proof  of  Theorem  (2.2),  Markov  property,  to  the 

process  p  - >Y(p,»),  which  has  exactly  the  same  form  as  the  process 

p  — >  Xp  of  (2.2) ,  we  see  that 


(4.28) 


yg  3  YGnH  [  YH 


Since  G  DG  HK  ,  (4.25)  and  (4.28)  yield 


(4‘29)  yg  =  yg  v  ygdk  1  ygah  v  ygdk  [  yh  v  yghk 


By  the  independence  of  the  increments  of  M  over  disjoint  sets,  we  have 


\  1 1  yq  v  yh  • 


which  yields,  when  used  with  (4.26)  and  (4.29), 


(4.30)  YG  ]  yGnH  v  YgDK  1  YH  V  YGOK  V  MK  =  YH  V  YK 


Note  that  the  values  of  Y  on  G  U  H  U  K  are  determined  by  M  on  that 
set,  whereas  the  values  of  X  on  L  are  determined  by  M  on  L.  So, 


XL  > 1  YG  V  YH  '  YK  ' 


which  together  with  (4.30)  and  (4.26)  implies  that 


(4.31) 


YlY  V  Y  TYVYVX. 

*G  J  GDH  GDK  1  H  K  L 


Since  X  =  X  v  X„_„,  we  may  put  X_,.  _  into  every  term  in  (4.31)  by  using 

t  jj  g  n  l  n  t 


(4.25),  which  gives 


(4.32) 


V  X 


ghl 


Y  V 

GHH 


GflK 


v  X 


GHL 


(  Y„  V  Y„  V  XT 


It  follows  from  (4.27)  defining  Y  and  (4.14)  for  X  that  X(p,t)  =  e_bt  X(p,0 
+  Y(p,t)  .  Thus,  X  on  GHL  together  with  Y  on  G  determine  X  on  G,  that  is. 


(4.33) 


XGHL  =  XG 


Putting  (4.33)  into  (4.32),  noting  that  XG  =  X^V  and  applying  (4.25) 

once  again,  we  obtain 


(4'34>  XG  1  YG nH  '  'WK  '  XGHL  '  XG n K  IYH  V  YK  V  *L  '  XGOK 


Finally,  noting  that 


G()H 


v  Y 


GOK 


v  X 


GHL 


V  X 


Gf)K 


=  X 


ghl 


XGHK  V 


XGHK  X9G 


t 

I 


and  that 


YH  V  YK  V  XL  V  XGHK  "  *L  V  XH  V  X  K 


we  see  that  (4.34)  is  in  fact 

XG  ]  X9G  [  XLUHUK  ' 


which  is  the  desired  statement  of  the  Markov  property. 
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